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ABSTRACT 

The observed apparent velocities and luminosities of the relativistic jets in 
AGNs are significantly different from their intrinsic values due to strong special 
relativistic effects. We adopt the maximum likelihood method to determine si- 
multaneously the intrinsic luminosity function and the Lorentz factor distribution 
of a sample of AGNs. The values of the best estimated parameters are consistent 
with the previous results, but with much better accuracy. In previous study, it 
was assumed that the shape of the observed luminosity function of Fanaroff- Riley 
type II radio galaxies is the same with the intrinsic luminosity function of radio 
loud quasars. Our results prove the validity of this assumption. We also find that 
low and high redshift groups divided by z = 0.1 are likely to be from different 
parent populations. 

Subject headings: galaxies: active — galaxies: luminosity function — galaxies: 
jets — galaxies: statistics 



1. Introduction 

The large scale and relativistic jets are important characteristics of active galactic nuclei 
(AGNs). Determining their proper velocity distribution and intrinsic luminosity function 
is of fundamental importance for understanding the central engine of AGNs. However, the 
observed apparent velocities and luminosities of the relativistic jets in AGNs are significantly 
different from their intrinsic values due to strong special relativistic effects, i.e. (3 app = fzf^g 
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and S = J^I^ (1 + z) 1+a , where @, (3 app , 9, L, S, 5 and cLl are the proper velocity, 
the apparent velocity, the angle between the axis of jet and the line of sight, the intrinsic 
luminosity, the observed flux, the Doppler factor and the luminosity distance, respectively. 
The index p is equal to n — a , where a is the spectral index (S oc v a ) and n = 2, 3 for 
continuous or discrete jet, respectively. 

Due to the lack of reliable and accurate methods to obtain (3 or 9 of each source, it 
remains difficult to determine directly the intrinsic luminosity and the proper velocity of a 
relativistic jet. Several authors have addressed these problems with statistical methods. For 
example, Padovani & Urry (1992) determined the evolution parameters of different types 
of galaxies with the < V/V m > method, and then used the observed luminosity function of 
Fanaroff- Riley type II radio galaxies (FR II) as the intrinsic luminosity function of relativistic 
jets in AGNs, after multiplying a constant factor. This approach is at best an approximation 
only, because the observed luminosity function of FR II may or may not resemble the intrin- 
sic luminosity function of relativistic jets in AGNs. Lister & Marscher (1997) extended this 
work with a simulation approach by fully considering the flux limit and Doppler beaming 
effect. However, they utilized the parameters determined by Padovani & Urry (1992) due 
to the size limit of the simulation. Therefore, the above mentioned approaches could not 
determine simultaneously the intrinsic luminosity and Lorentz factor distribution of rela- 
tivistic jets in AGNs. Actually, the observed luminosity function and the apparent velocity 
distribution can be reproduced if the intrinsic luminosity and the Lorentz factor distribution 
are given (Lister 2003, Vermeulen & Cohen 1994). Thus the intrinsic luminosity and Lorentz 
factor distribution could be inferred simultaneously with the maximum likelihood method. 
Therefore, we could test the previous assumption that the observed luminosity function of 
FR II has the same shape with the intrinsic luminosity function of the radio loud quasars. 

The form of the maximum likelihood estimator could be written as 



where % denotes the number of each source in the sample and p is the probability density of 
detecting a source for the given apparent velocity, observed luminosity and redshift. C is 
the normalized factor. The probability density p is 



where /3 apP , d-A, T(z), f((3 app \ L Q , z) and &(L ,z) are the apparent velocity, the angular 
distance, the look-back time, the apparent velocity distribution (for the given observed lu- 
minosity and redshift) and the observed luminosity function, respectively. As proved by 
Vermeulen & Cohen (1994, also see Appendix A at the end of this paper), if we assume 



MLE = In p(P apP ,i, L a , 



Zi), 



p(p app ,L ,z) = f(/3 app \ L , z)$(L , z)d A (zf(l + zfc 
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the intrinsic luminosity function is a single power law (<&(L) oc L~ A ), L Q will not appear 
in f(Pa PP \ L , z). Therefore, we could deal with the apparent velocity distribution and the 
observed luminosity function separately 

We present the results obtained from the apparent velocity and the observed luminosity 
(also the redshift) in §2. The implication of our results for the unification scheme of radio 
loud AGNs is addressed in §3. In §4 we summarize our conclusions and make some discus- 
sions. The cosmological parameters adopted throughout this paper are (H , Qm, ^a) = 
(75 km/s/Mpc, 0.27, 0.73) . 

2. Sample selection and data analysis 

To investigate the real parent population of the radio loud AGNs, we should adopt a 
large sample containing different kinds of radio loud AGNs, including radio galaxies, quasars, 
and blazars. The largest available sample of the kinetics of relativistic jets in AGNs is from 
the "VLBA 2 cm survey" (Kellermann et al. 2004). Its successor, the "MOJAVE survey 0", 
provides a complete flux-limited sample further. Due to the high frequency observation (15 
GHz), any radio emission from large scale structures is effectively excluded. The majority 
of the sample are the flat-spectrum radio quasars (FSRQ), which is the most luminous type 
in AGNs. According to the unified model, FSRQ is the beamed version of FRII due to the 
small angle between the jet and the line of sight. 

2.1. Analysis of the apparent velocity data 

In this section we only make use of the apparent velocity data to determine the intrinsic 
luminosity and Lorentz factor distribution. The detailed formulae of the apparent velocity 
distribution are presented in Appendix A. 

Our sample of the apparent velocity data comes from the "VLBA 2 cm survey" (Keller- 
mann et al. 2004). We only adopt the fastest component of each source that has a quality 
factor of "Good" or "Excellent" (see Figure 11 in Kellermann et al. 2004), since several 
authors have pointed out that the velocities of some patterns are much slower than the bulk 
velocity responsible for the Doppler boosting effect (Cohen et al. 2006). The sample used 
here contains 16 BL Lac objects, 12 radio galaxies and 76 quasars. We assume the Lorentz 
factor distribution is also a single power law (^(7) oc 7° (1 < 7 < 7 max )) as adopted by 
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previous studies. As indicated by previous studies, the maximum value of f3 app is nearly the 
same as 7 maa; (Vermeulen & Cohen 1994). The maximum value of f3 app in the sample is about 
31, and therefore we set 7 max = 32 for safety. As assumed by the previous study (Lister & 
Marscher 1997), we set n = 2 and a — 0. It has been found that these simple assumptions 
could give reasonable results. The estimated result from the whole sample is shown in Figure 
1. Obviously, the constraint on the parameters is weak, which is due to the limitation of 
the observed data. However, it is still useful to constrain the values of parameters with such 
data, because this allows us to compare between the results from the apparent velocity data 
and that obtained from the observed luminosity data, in order to check the consistency of 
our models. 

By fixing the index A, we could investigate whether the Lorentz factor distribution 
evolves with the redshift. We fix A = 3.10 (the best estimated value of A for the whole 
sample) and estimate a in each redshift bin. As shown in Figure 2, the Lorentz factor 
distribution is much steeper in the first bin (0 < z < 0.1) than in other bins. The difference 
of the apparent velocity distributions between the first bin and other bins is also significant 
at the 99.98% level obtained by K-S test. Therefore, we should divide the sample into low 
and high redshift groups by z = 0.1. The estimated parameters of the two groups are also 
shown in Figure 1. The difference between the results from the whole sample and the high 
redshift sample is only slight, whereas the confidence region of low redshift sources differs 
with others at 90% confidence. 



2.2. Analysis of the observed luminosity function and redshift distribution 

We assume the local intrinsic luminosity function is a single power law 

$(L, z = 0) = kL' A (Li < L < L 2 ), 

and to compare with the previous results, we also assume a pure luminosity evolution model 
as used in Padovani & Urry (1992) and Lister & Marscher (1997), 

z) = $(L/e X p[T(z)/T], z = 0), 

where r is the evolution parameter in unit of Hubble time. 

Using the equations of the beamed luminosity function (see details in Appendix B), 
we can estimate the parameters in the intrinsic luminosity function, the Lorentz factor 
distribution and the evolution form by the maximum likelihood method. We combine the 
data from "VLBA 2 cm survey" and "MOJAVE survey" (Figure 3). We again divide the 
sample into low and high redshift groups by z = 0.1. The values of the best estimated 
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parameters from high redshift group are a = — 1.73±S;||, A = 2A8tom, r = 0.29±°;oi and 
log(Lx(W/Hz)) = 26.25l°;°2; the results from the low redshift group are a = —1.5+^, 
A = 2.8ti° and log(Li(W/Hz)) = 24.39±g;2o ( no evolution is assumed for the low redshift 
group, and la error for each parameter of concern is shown). The value of log(L 2 (W/Hz)) 
is fixed at 34, since higher values barely change the results. Figure 4 shows the confidence 
region for (a, A) by fixing other parameters at the best estimated values. Although the 
confidence regions in Figure 4 for the low and high redshift groups are consistent within the 
68% confidence region, the lower intrinsic luminosity limits of the two groups are different 
significantly. As also shown in Figure 4, the result from the luminosity data for the low 
redshift group is only marginally consistent with the result from the apparent velocity data 
within the 68% confidence region. This indicates that our current model do not describe 
the low redshift group very appropriately. Actually, if we adopt the whole sample, the 
result of the Lorentz factor distribution is unreasonable, i.e. a > is required to obtain 
the maximum value of the estimator, which is ruled out by previous studies. For example, 
if we adopt log(Li(W/Hz)) = 24.39 indicated by the low redshift group, the values of best 
estimated parameters of the whole sample are a = 30.6, A = 3.5 and r = 1.0. The K- 
S test is performed to estimate the goodness of fit; the corresponding probabilities of the 
apparent velocity, the observed luminosity and the redshift distribution are 10~ 22 , 10% and 
7%, respectively. Therefore, the consistency is quite poor, especially for the apparent velocity. 
If the lower intrinsic luminosity limit is adjusted as a free parameter, the value of a will 
become larger and even without an upper limit, which is more inconsistent with the result 
from the apparent velocity data. Therefore, we find that the results from both the apparent 
velocity and the observed luminosity data indicate the low and high redshift groups are likely 
to be from different parent populations. Due to the limited size of the low redshift group, 
we only utilize the high redshift groups in further analysis. 

As shown in Figure 5, there is slight difference between the result obtained from the 
apparent velocity data and the observed luminosity function (within 68% confidence region). 
This is mainly due to the simple evolution form used here. When performing the K-S test 
to estimate the goodness of fit, we find the result of the best estimated value of parameters 
could only be accepted marginally. The corresponding probabilities of the apparent velocity, 
the observed luminosity and the redshift distribution are 45%, 56% and 6%, respectively. 
However, if we adjust the value of r to the la boundary (r = 0.32), the result of K-S test 
could be improved. The corresponding probabilities of the apparent velocity, the observed 
luminosity and the redshift distribution are 45%, 61% and 20%, respectively. The effect of 
a higher flux limit is only mild (see the inset in Figure 4). 
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3. Implication for the unification scheme of radio loud AGNs 

There are two key parameters in the unification scheme of radio loud AGNs. One is the 
orientation of the sources, and another is the intrinsic luminosity of the sources. The high 
luminosity population contains quasars and luminous radio galaxies (FR II), while the low 
luminosity population contains BL Lac objects and less luminous radio galaxies (FR I). On 
the other hand, the viewing angles of quasars and BL Lac objects are smaller than that of 
radio galaxies. Therefore, the superluminal motion is more common in the aligned objects, 
and the observed luminosity is strongly affected by special relativistic effects. However, the 
relativistic beaming effect of radio galaxies is mild and the observed luminosity function of 
radio galaxies is supposed to be similar to the intrinsic one. The discussion about the low 
luminosity population could be found in Padovaani & Urry (1990, 1991) and Urry et al. 

(1991) , and the unification scheme was reviewed in Urry & Padovani (1995). Due to the 
limited sample size of low luminosity population and unclear relationship between FR I and 
FR II, we focus our discussion on the high luminosity population below. 

Several papers investigated the luminosity function of AGNs. For example, Urry & 
Padovani (1995) found the evolution parameter r of FR II and FSRQ are 0.26+q; 10 and 
0.23^q;q4, respectively. Padovani & Urry (1992) adopted a double power law to fit the 
observed luminosity function of FR II, which was identified as the parent luminosity function. 
The power law indices of low and high luminosity band were found to be B 1 = 2.48+°;}^ and 
B2 = 3. 9lg 7> respectively. Here we have determined simultaneously the intrinsic luminosity 
function and the Lorentz factor distribution of relativistic jets in AGNs by the maximum 
likelihood method. However, we find a single power law form of the intrinsic luminosity 
function is sufficient to describe the sample we used here. For comparison, we have re- 
analyzed the observed luminosity function of the FR II sample adopted in Padovani & Urry 

(1992) by the maximum likelihood method (for completeness, only the sources with z < 1.3 
are used). To compare with our result, we also adopt a single power law ($(£,,, z = 0) oc L~ A ) 
and the pure luminosity evolution. The results are A = 2.5llo!i5 an d r = 0.40lo!o9- The 
corresponding probabilities obtained by K-S test of the observed luminosity and the redshift 
distribution are 75% and 80%, respectively. Therefore, we find this simple model could 
describe the data well, and the result is consistent with the intrinsic luminosity function we 
have obtained in §2. In previous studies (Padovani & Urry 1992; Lister & Marscher 1997), 
it was assumed that the shape of the observed luminosity function of FR II is the same with 
the intrinsic luminosity function of the radio loud quasars. Therefore, our results prove the 
validity of this assumption. The single power law form of the luminosity function may be 
somewhat simplified, and a more complex form will be investigated when a larger and more 
complete sample is available. 
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4. Discussions and conclusions 

We have determined simultaneously the intrinsic luminosity function and the Lorentz 
factor distribution of relativistic jets in AGNs by the maximum likelihood method, and have 
confirmed the previous assumption about the shape of the intrinsic luminosity fucntion. 
The result of the Lorentz factor distribution is also consistent with the previous result. For 
example, Lister (1997) claimed the index of the Lorentz factor distribution a is roughly 
between -1.75 and -1.5. However, he fixed the parameters of the intrinsic luminosity as a 
prior condition. 

In the previous studies, it has been assumed the Lorentz factor distribution is the same 
at all redshift. However, we find the Lorentz factor distribution is much steeper at low 
redshift (z < 0.1) with the result from apparent velocity data, though the uncertainties of 
the results are large. This indicates that the low and high redshift groups are likely to be 
from different parent populations, i.e. the dual-population scheme (Jackson & Wall 1999). 
The majority of low reshift sources are low luminosity ones. They are not as energetic as the 
high luminosity quasars. Therefore, the most extremely relativistic jets are relatively rare in 
this population. However, this is not indicated by the results from the observed luminosity 
function. As discussed in §2, this is likely to be due to the relatively simple model applied 
here. Further more detailed study could be performed when the sample is large and complete 
enough. 

We assume the Lorentz factor is independent of the intrinsic luminosity, i.e. the 
L — 7-independent (LGI) model. Lister & Marscher (1997) investigated a particular L — 7- 
dependent (LGD) scenario (L oc 7^). They found the predictions of the best-fit LGD model 
were very similar to the best-fit LGI model but predicted very few high viewing angle sources 
compared with the Caltech-Jodrell Bank sample. We speculate this may be due to the some- 
what arbitrary form of the relation between the Lorentz factor and the intrinsic luminosity. 
Actually, even when the LGI model is employed, there should be some correlation of the 
Lorentz factor and the intrinsic luminosity in the observed sources due to the higher lumi- 
nosity threshold at higher redshift. 

Besides the shape of the intrinsic luminosity function, the space density of different 
types of AGNs should match with each other according to the unified scheme. However, 
both of the space densities of FR II and FSRQ evolve with redshift significantly, and the 
value of the dividing angle of different types of AGNs may also be related with the luminosity 
and redshift (e.g. Willott et al. 2000; Arshakian 2005). It is therefore quite complicated to 
demonstrate the consistency of space densities of different types of AGNs within the context 
of the unification scheme, and thus beyond the scope of the current work. 
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A. Calculation about the apparent velocity distribution 

The approach presented here is similar to that in Vermeulen & Cohen (1994). The 
apparent velocity probability density function (pdf) p(j3 app ) could be obtained by the differ- 
entiation of the cumulative distribution function P(f3' app < Papp), 

p(Pa PP ) = dPiP Z ~ ^ , (Al) 

Q'Papp 



and 

P(P' apP < Aw) = 1 - P(P' app > Aw) = 1 - / dry I p(9, 7 )d0 , (A2) 



where 

ll = \/Pl PP + 1 , 

1u Tmax i 



cose! = /3 ° pp7+v/ °' 2 ii £ , 

f 0, p apP <P, 

COS 6 U = I /3g pp 7-y/7 2 -l-/3g pp R R 

( (1+^)^ ' Pa ™-^ 

and 7 = l/y/l -p 2 . 

Due to the Doppler boosting effect, p(9, 7) depends on the observed luminosity L Q 
LS n . We assume pdf of L has a power law form, 

p(L) = 



(-^mfn ^max) 
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then 

p(L \ 6 }1 )dL = p(L = [7(1-/3 cos 6)} n L )dL 

_ q{[-y(l-l3cose)] n L }-(i+ 1 ) 



, , „ -[ 7 (l-/3cos0)] n dL o 

^min ^max; 

g[ 7 (l-/?CQBfl)]-"g ,r 



Therefore, 

p tt/2 



p(0, 7 | £„) =p(L | ^7M%(7)/[/ 1 7ma >(7)^7/ c ^ p(0)p(L o |0, 7 )d0] 

= [7(1 - /? cos g)]-^ sin ^( 7 )/[/ 1 7max mt)^ [7(l - /3)r ;;_ + r 7 

Note that there is not L Q in equation (A3). Therefore, 

p(0, 7 |L o )=p(0, 7 ). 
Substitute equation (A3) into equation (A2), we have 



> Aw) = J 7 ^7 l)d6 

= C/ 7 ]>(7)d7£° s t[7(l-^)]-^ 



(A3) 



= ^ p« p(t) rf 7(1-/3 cos gp] - [7(1-/3 cos 

^7! 7/3 ' rig— 1 

= C/J* {^[7(1 -/3cos^)]- n9+1 - [7(1 -/9cos^)]^ +1 }rf 7 • 

To be consistent with the notation in Appendix B, we denote A — q + 1. 

Using equation (Al), we could obtain the pdf of the apparent velocity distribution. In 
this paper we assume ^(7) is a power law, i.e. p(j) = C^ a . The examples of different values 
of a are shown in Figure 6 



B. Calculation about the apparent luminosity function 

The method presented here is the same with that in Lister (2003). Here we still list the 
main equations for clarity. 

To calculate the apparent luminosity function, we should know pdf of 5. We assume the 
orientation of jets is random, i.e. p{9) oc sin(0). Therefore, we have 

Jc(8) V7 2 - 1 

where 

8-\ 1-L<S<1, 

1< 5 < 5 = 'y (3 + 1) 

2£ ) ± — u — "max /maxVPmax ~ ^-J 5 



c(S) 
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We assume the form of the intrinsic luminosity function is 

<f)(L) = kL- A (Li < L < L 2 ) . 

Since 

$(L o ,0) = J P (L ,5)d5 = J p (5)<f>(L)^-d5 , 

and L Q = L5 n , we have 

$(L ) = kL- A / 5- n+nA - 2 d5 / -^J=d 1 , 

where 

= 7 m L &(A0 = (^o/^i) 1/n (£i7 m : x < A, < iiO » 
a(L ) = 7-^, b(L ) = 5 max (^iC ax < L < L 27 ^ x ) , 

a(L c ) = (L C /L 2 )V-, 6 (l o ) = £ max (L a7 - < L < L 2 < ax ) , 

In this paper we assume p( 7 ) is a power law, i.e. p(j) = C 7 a . The results of different 
values of a are shown in Figure 7. This figure is similar to the Figure 3 in Lister (2003) (We 
mention in passing that there are some typos in the figure caption.). 
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Fig. 1. — Inferred parameters from the apparent velocity data. The confidence regions are 
results obtained with the whole sample (solid lines), z > 0.1 sources (dashed lines), and 
z < 0.1 sources (dot lines), respectively. The 68% and 90% confidence regions are shown for 
each case. 
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Fig. 2. — Evolution of the index of the Lorentz factor distribution by fixing A = 3.10. From 
low to high reshift, there are 13, 17, 17, 13, 17, 17 and 10 sources in each redshift bin, 
respectively. The la errors are shown for each redshfit bin. 
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Fig. 3. — The luminosity and redshift sample (MOJAVE survey). The solid and dashed lines 
are the flux limit corresponding to 0.32 Jy and 0.50 Jy, respectively. For completeness, only 
the sources with z < 2.5 are included. 
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Fig. 4. — The results from the observed luminosity and redshift (solid lines). For comparison, 
the results from the apparent velocity are also shown (dashed lines). The 68% and 90% 
confidence regions are shown for each case. 
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Fig. 5. — The thin solid lines are the results from the observed luminosity and the redshift 
distribution. The dashed lines are the results from the apparent velocity distribution (z > 
0.1). The thick solid lines are the combined results. The inset shows the results from the 
observed luminosity and the redshift distribution for S min = 0.32 (solid lines) and S min = 0.50 
(dashed lines). The 68% and 90% confidence regions are shown for each case. 
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Fig. 6. — The pdf of the apparent velocity distribution with ^(7) oc 7 a (l < 7 < 30) and 
n — 2. Three cases are shown, i.e. a = — 1 (solid line), a = — 2 (dashed line) and a = — 3 
(dot line). 
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Fig. 7. — The observed luminosity function with p("f) oc 7 a (l < 7 < 30) and n = 2. The 
thick solid line is the intrinsic luminosity function, i.e. 4>(L) oc L -2 (Li < L < 10 7 Li). Three 
cases are shown, i.e. a = — 1 (thin solid line), a = — 2 (dashed line) and a = — 3 (dot line). 



